On the "effective exponent theory" of the Coulomb Luttinger liquid 
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The "effective exponent theory", developed by Wang, Millis and Das Sarma in [Phys. Rev. B 69, 
167101 (2004), Phys. Rev. B 64, 193307 (2001)], fails to calculate correctly the dynamic correlators 
of Coulomb Luttinger liquid. Main drawbacks are (i) cutting off the Coulomb potential by killing 
its long-range component, (ii) the absence of non-contradictory procedure to determine the 'scaling 
cutoff' and correlation functions for true Coulomb interaction. 



I. INTRODUCTION 



II. COULOMB INTERACTION PROBLEM 
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Wang, Millis and Das Sarma (WMS) proposed^ an 'ef- 
fective exponent' approach to calculate the dynamic cor- 
relation functions in a Luttinger liquid with Coulomb in- 
teraction (CLL) and developed this theory in recent pub- 
lications^ opposing it to our approach^*& to the same 
problem. They claimed that (i) their method is 'more ele- 
mentary' than our one, but is applicable in a wider energy 
range, (ii) both methods are equivalent in describing the 
threshold behavior of charge-density wave (CDW) struc- 
ture factor and the dynamic spectral function (after a 
technical error in deriving the spectral function in our 
work4 has been eliminated^). 

In this Comment, we show that: 

i) The 'effective exponent theory' fails to take ad- 
equately into account the long-range nature of the 
Coulomb interaction and does not correctly describe the 
dynamic correlation functions. 

ii) WMS's method contains a concealed procedure of 
fitting to the results of independent approaches and hence 
is not a self-sufficient theory. Correspondingly, the range 
of applicability of WMS's method is at least not wider 
than that of the reference approach. 

iii) WMS's results^ for the CDW structure factor and 
spectral function disagree with our asymptotically ex- 
act results^i^, and also contradict their own earlier nu- 
merical calculations^. In particular, the ratio of two re- 
sults for the spectral function p(q, to) tends to infinity 



in the threshold region of 8 

13/2 



-0 like 



exp[C| lni5| 3 / 2 ], that is faster than any power of 1/5. 

The theory of WMS would not probably have been 
deserving a special discussion, if it had not involved an 
important physical problem of Coulomb interaction ef- 
fects in low-dimensional electron systems and methods 
of its description. 

Coulomb interaction effects are a rather complicated 
matter, studying which requires a delicate approach. The 
brash approximations can destroy the studied model so 
that it will be far from the realistic physical system. The 
analysis of the WMS theoryi2^ is quite instructive to il- 
lustrate this. Their theory distorts strongly the Coulomb 
potential by killing its long-range component that is of 
most interest. 



In principle, the Coulomb interaction problem is solved 
in the spinless Luttinger liquid theory. The bosonization 
method allows one to exactly calculate all the correlation 
functions for any type of the electron-electron interaction 
potential, and, in particular, does so for the long-ranged 
Coulomb interaction!^ As a result, correlators are ex- 
pressed in terms of definite integrals. Thus, the CDW 
structure factor <Scdw(<Zi ^) equals^iifl 



S, 



CDW 



dx e 1 



dte luJt S(x,t) , (1) 



with 
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-^2 d l ( exp[-u F / — 



+ OC 



dp 



(1 _ e -^ P t-^) e -«|p|] C0S (2fc F x)) . (2) 



Here vp is the Fermi velocity, oj q = \q\vp / g(o) is bosonic 
excitation spectrum, the interaction parameter g(q) is 
related to the Fourier-component of the interaction po- 
tential by g(q) = [1 + V(q)/hTrv F ]~ 1/2 . 

Therefore, the task is reduced to the calculation of such 
integrals. However, in the case of Coulomb interaction 
such calculation entails serious mathematical difficulties 
that are connected inherently to the slow decrease of the 
interaction potential at large distances. 

True Coulomb interaction is well known to be a special 
one. The distinctive feature of Coulomb interaction is 
that the integral of Coulomb potential 



V(x) 



d 2 



diverges at large distance, d being the quantum wire di- 
ameter. Correspondingly, the Fourier-component of the 
potential goes to infinity in the long-wave limit q — > as 



V{q) ~ 2e 2 \lnqd\ 



(3) 



This manifests itself in the logarithmic divergency of the 
collective excitation velocity v(q) at q — > 0. Any approxi- 
mation of the Coulomb interaction by cutting off at large 
distance loses this feature and hence makes the interac- 
tion a non-Coulomb one. 
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The central question is whether long-ranged Coulomb 
interaction leads to a qualitatively new behavior of ID 
electrons as compared to the short-ranged interaction 
case, which is easier to solve? 

The classic paper of H.J. SchulaS gave the positive an- 
swer to this question. Coulomb interaction was shown to 
have dramatic consequences for the density-density cor- 
relations of a ID electron liquid. Specifically, the pair 
density correlator in the Coulomb case decays with the 
distance in the extremely slow fashion exp[— c(lnx) 1 / 2 ], 
instead of a power law, typical of the short-ranged in- 
teraction. In other words, in ID systems Coulomb in- 
teraction induces a quasi-long-range order, close to the 
Wigner crystal. 

Ref.0 deals only with the static correlations. Dynamic 
correlations are much more difficult matter because the 
dependence on both q and u> has to be addressed. Exact 
results are apparently of key importance in the explo- 
ration of this problem. They were not obtained until 
very recently. 

The analytic approach to dynamic CLL correlators de- 
veloped in our paper— consists in the regular derivation of 
the integral equation for ScuwiljUj) and p(q,uj) and its 
solution near the threshold. What is important is that 
no model simplifications pertaining to the form of the 
Coulomb potential are used in this approach. 

In WMS's method^* 2 * 3 ., no Coulomb integrals are cal- 
culated analytically. WMS replaced the Coulomb inter- 
action with an 'effective' short-ranged interaction, i.e. 
cut off the long-ranged Coulomb potential at large dis- 
tance. Thus their approach is a model one. But this is 
not its major shortcoming. The real trouble lies in the 
arbitrariness of such replacement, as a consequence of 
which the authors, turning to their method in different 
publications^* 2 *^, each time obtain different results. This 
becomes evident if one considers how Coulomb interac- 
tion is transformed into short-ranged one. Since WMS 
in their publications do not go into the details of such 
procedure, we briefly outline it here for the case of the 
CDW structure factor Scdw{q,uj). 



III. EFFECTIVE EXPONENT THEORY 

WMS starts with imposing the long-distance cutoff on 
the Coulomb potential at some wave- vector go, 



teraction, according to which 



V(q) = 2e 2 \H(q Q +q)d)\ 



(4) 



In contrast to Eq. ©, this potential does not diverge 
at q — * 0, and hence is not the long-ranged Coulomb 
potential. 

The structure factor S'cdw(<Z, uj), having been calcu- 
lated with the screened potential of Eq. |0J, acquires the 
dependence on the cutoff parameter, i.e. Scuwiq, uj) is 
transformed into 5cdw(9o; 3,w). However, even for such 
cut off interaction no formulas are derived. Instead, the 
authors refer to common results for the short-ranged in- 



Stt)w(q,u) oc (uj - LL> q -2k F ) ac 



(5) 



where the exponent equals ocdw = ff(0) — 1. 

WMS 'simply replace' the Coulomb interaction param- 
eter g(q) ~ | In |gd|| -1 / 2 by g(qa) and postulate that the 
'effective' exponent equals 



acDw(go) = g(qa) - i 



(6) 



Determining go is the main trouble of WMS's ap- 
proach. Indeed, no consistent theory exists until the pro- 
cedure to determine go is defined. Assuming go constant 
(a fitting parameter) results in a standard power-law de- 
pendence of the structure factor Scbw(<1,uj). For CLL, 
this corresponds neither to the numerical calculations nor 
the expectations, based on exactly solvable models. To 
obtain the formulas similar to those of Schulz, 8 WMS 
subject 5cDw(go; 9i w ) to further modification by intro- 
ducing the explicit dependence of go 011 uj an d g. In other 
words, the cutoff parameter go is replaced by a function 
of two arguments qo(q,uj). Notice that introducing the 
cutoff go, which depends on uj and g, can not be justified 
by any renormalization group arguments. 

WMS propose no regular way to obtain the function 
go(g, w). Consequently, their method does not solve the 
Coulomb interaction problem. In fact, go is chosen from 
the intuitive conjectures, differently for different correla- 
tors being calculated. In the case of «ScDw(g, uj), WMS 
postulate without convincing justification that qo(q,uj) is 
determined from the following equation, 



UJq-2k f 



(7) 



In the case of p(q,uj), they state that go(g, uj) should be 
determined from other equation, 



UJ U 



Uq-k B 



(8) 



In doing so, they ignore the fact that in both cases 
go (g, uj) enters into the same expression for the potential 
in Eq. |0J. As a result, the electron-electron interaction 
potential turns out to be different for different correla- 
tors. This 'scaling cut-off' procedur e) 1 '^ 3 is seen to be 
physically inappropriate. 

The most important problem WMS confront further 
is how to obtain the end result for correlation functions, 
even having the function qo(q,u). In an earlier version 
of their theory^, <S*cdw was determined according to Eqs. 
©J©,©, but this result did not satisfy the authors, and 
in the final release 3 - they declared another procedure. 

WMS postulated that the structure factor satisfies the 
following differential equation (see Eq. (14) of Ref. 0) 



d In S, 



CDW 



(0 



91ne 



acD\v(<7o) 



(9) 



where e denotes a normalized energy deviation from the 
threshold, e = (uj — uj q ^2k F )/uJo, with ujq = vf/ fid, /3 = 
\KVFl2e 2 ] 1 ' 2 . 
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Eq. 10 is seen to be a fundamental principle of the 
'effective exponent theory', which allows WMS to find 
the dynamic structure factor analytically. However, no 
proof is provided that the structure factor, given by Eqs. 
lJU,©, indeed satisfies this equation, at least asymptot- 
ically in the threshold region of e — * 0. 

Solving Eq. @ , WMS obtain the following expression 
for Scow (Eq. (15) of Ref. H, 



5 C Dw(g,w) oc exp(-2/3|lne| 1 / 2 )/e 



(10) 



The asymptotically exact expression for 5cdw (<7i w ) & t 
e — > was found in our paper—, 



Scuwiq, w) 



vp e 



-4/3| lne| 1/2 



ellnelVa 



(11) 



The comparison of the WMS result (Eq. JT^Jl) with our 
result (Eq. iJTTjl') via dividing one expression by another 



WMS 

Sgs 



oc exp(2/3|lne| 1 / 2 ) 



(12) 



shows that in the most interesting threshold region e — > 
+0, the ratio goes to infinity. This means that above re- 
sults for the CDW structure factor are not equivalent, in 
contrast to what WMS state'. Possibly, one may invent 
some prescript to fit the effective exponent expression to 
the correct result but this hardly will be of interest. 

Exactly the same steps are used in 2 ^ to obtain the 
dynamic spectral function. The latest presentation of the 
effective exponent theory^ does not contain the explicit 
form of the spectral function. However, WMS's strategy, 
as it is presented in Ref. y, straightforwardly leads to the 
following expression, 



p^cxexph^lln^/ 2 ] 



(13) 



which gives the threshold behavior of the CLL dynamic 
spectral function in the spinless case. Our result for the 
spectral function^ is derived below, see Eq. (|19|) . In the 
threshold region S — ► +0, the ratio 



cx exp — m<y p 

PGS Off 



(14) 



goes to infinity faster than any power of 1 / 8, which means 
that the results are incompatible. 

It is interesting to note that the results of analytic 
effective exponent theory contradict WMS's earlier nu- 
merical calculations also in the spinful casei. Indeed, the 
analytic 'effective exponent procedure' of Ref. gives in 
this case p(q,u>) oc exp[— ln<5| 3 / 2 ], whereas in Ref.Q] 

it was found that p(q,u>) oc exp[— ^ | In t5| 3 / 2 ] . The ratio 
of these expressions again goes to infinity in the threshold 
region as 6 — > +0. 

The above examples show convincingly that the 'ef- 
fective exponent theory' is a fitting to the independent 
(analytic or numerical) results, rather than a justified 



self-sufficient procedure. Consequently, the WMS claim 3 
that their 'effective exponent theory' applies over a much 
wider energy range than our method is misleading. The 
effective exponent theory can not give correct results out 
of the applicability range of independent reference theory. 



IV. ON THE FITTING TO NUMERICAL 
CALCULATIONS 

As an independent reference theory in their earlier 
worki, WMS used the numerical calculation in the frame 
of standard Luttinger model. WMS chose the functional 
dependence of the cutoff parameter qo and correlation 
functions by fitting to the numerical results. For this pur- 
pose, the standard expressions for the correlation func- 
tions of the Luttinger model with short-ranged interac- 
tion were modified by introducing fitting parameters in 
the exponents. 

The shortcomings of such approach are evident. First, 
correlation functions are of primary interest in the thresh- 
old regions, where they either diverge or quickly go to 
zero; however, it is in these regions that the accuracy of 
numerical calculations is least trustworthy. Second, the 
fitting formulas do not follow from any physical theory, 
but rather represent a formal mathematical interpola- 
tion. 

Moreover, one has serious grounds to doubt whether 
WMS's fitting to numerical calculations is valid. 

1) The formula for the electron Green function used 
in WMS's calculations is incorrect. Specifically, Eqs. 
(4), (5) of Ref. contradict to the very definition of 
the Green function. Green's function is defined ini as 
G r (x,t) = (ip r (x, t)ipj,(0, 0)), where r refers to the right- 
moving (r = +1) or left-moving (r = —1) particles. The 
Fourier-transform of the Green function G r {q,uj) pos- 
sesses the fundamental property: at zero temperature, 
G r (q,uj) = if uj < 0, the energy u> being measured 
from the chemical potential. 1 - WMS's results (4), (5) 
show that the Green function of the left-moving par- 
ticles G T —-i(x, t) is analytic in the upper half-plane of 
the complex time- variable. This immediately gives that 
Green's function G r= _i(g, w) ^ at id < 0, but instead 
G r= _i(g, ui) = at uj > 0. Thus analytic and spectral 
properties of the electron Green function are violated in 
WMS's worki. 

2) Fitting procedure uses the behavior of the spec- 
tral function p(q, uj) at frequencies much larger then the 
Fermi energy, where the Luttinger model is not justified. 
Thus, they introduce some 'previously overlooked energy 
scale' oj s that is defined as uj s — 50(vF/d)y/Vo exp(l/Vb), 
with Vq = 4e 2 /irvF being a dimcnsionless measure of the 
interaction strength. The spectral function is considered 
at u > uj s . Note that uj s is of several orders of magni- 
tude greater than the electron Fermi energy, which is far 
beyond the scope of the exploited approach! 

3) The expression for the quasi-lD Coulomb interac- 
tion potential is incorrect in the high-momentum region, 



4 



where analytical results were fitted to numerical ones. 
WMS took 



V(q) = 2e 2 In 



A 



(15) 



with the short-range cutoff A w 2.5/ d, d being the diam- 
eter of the quantum wire, see Eq. (2) of Ref. U WMS 
denote the short-range cutoff also by qo, but here we re- 
serve qo for the long-range cutoff. WMS claim that the 
potential should be '1/q for q larger than some scale A set 
by the geometry and the wave function size'ii Eq. H15|) 
indeed guarantees that. 

However, such behavior is incorrect, since the potential 
should behave as 1/q 2 in the high- momentum region 12 . 
The correct formula is derived in Appendix A, 



V(q) 



In 



A 2 



(16) 



It describes adequately the behavior of the quasi- ID 
Coulomb potential both in the long-wave (V(q) « 
2e 2 |lng<i|) and in the short-range (V(q) ~ e 2 A 2 /q 2 ) re- 
gions. 

Generally, WMS uncritically overestimate their results 
pertaining to the region of high q and lu. 

Thus, using Eq. (|15|) . WMS found the strange behavior 
of the plasmon frequency Lu q at q ^> d , namely Lu q ~ 
qvp + const. This is not a unique property of CLL, as 
they assert, but a trivial mistake. Correct Eq. (|16fl gives 
Ljq ~ qvF, without any additional constant. 

They also discovered some dramatic modifications of 
p(q, w) at w > lu s as compared to the short-ranged inter- 
action case due to 'the slow (1/q) decay of the Coulomb 
interaction in the large momentum region'. The passion 
fades if one remembers that Eq. I|15|l for Coulomb poten- 
tial is incorrect, and bosonization does not work there. 



The conservation laws impose the following restriction 
on the spectral function: p(q,co) = at lu < u> q ; thus 
p(q, lu) has a threshold at u — > cu q + 0. 4 Applying this ar- 
gument to the RHS of Eq. {T7J|, we find that the integrand 
is nonzero in two narrow bands of the width S/dy/\ In S\, 
where S — (uj—uj q )/ujo, situated near to Q « and Q w q. 
The direct calculation shows that the contribution of the 
band near to Q « q to the integral is smaller than the 
one from the band Q « by the factor of S. There- 
fore, considering the asymptotic behavior of p(q,u>) at 
5 — > 0, we can retain only the region of Q « 0. On these 
grounds, we expand the RHS of Eq. fTTjl in the powers 
of Q, appearing in the combination (q — Q), 



p(q,w) 



4lu 



+ OC 



dQ 



1 



<?(<?) 



p(Q,u 



In the case of Coulomb interaction, it is sufficient to con- 
sider only the first term of the expansion in the RHS 
of Eq. i|18fl to obtain the correct asymptotic behavior of 
p(q,w) at £-> +0i£J£ 

On the RHS of Eq.JTg), the integral f+™ dQp(Q, lu ~ 
Lu q ) equals, up to a factor, the tunnelling density of states 
N(lu — u>q), which has a soft gap in the Coulomb inter- 
action case. Consequently, our method shows that the 
spectral function p(q,Lu) turns to zero as lu — > lu q , i.e. 
it contains a pseudogap. The functional dependence of 
the tunnelling density of states can be obtain either via 
direct integration, or, more elegantly, by solving the in- 
tegral equation, similar to (|17fl . 

Finally, in the limit S — > we obtain^ 

Alln^l 1 / 2 



p(q,u) 



.exp[-^|ln<5| 3 / 2 ] 



(19) 



where A(q) = [l — g~ 1 (q)sign(q)] . The obtained form 
of the spectral function is consistent with numerical and 
experimental results 1 ^* 1 ^ 1 ^* 1 ^ without any mystique. 



V. THE ANALYTIC APPROACH TO CLL 

In our papes^, we have developed an analytic method 
to determine the exact behavior of the dynamic corre- 
lation functions of the Coulomb Luttinger liquid in the 
vicinity of the threshold. We have taken into account 
that the velocity v(q) of bosonic excitations in a CLL di- 
verges as q — > and effectively used this fact. We will 
now briefly describe it to show how to obtain the correct 
expression for the spectral function on a regular basis. 

As it is shown in£, the r-fermion dynamic spectral func- 
tion satisfies the following integral equation 

v F f +cc 

p(q,u) = J dQK.(q-Q)p(Q,LU-u> q - Q ) , (17) 

where the kernel is /C(£) = [r — <7 -1 (£)sign(£)l 2 , the wave 
vector q being measured from the Fermi wave vector. To 
be specific, consider right-moving particles. 



VI. CONCLUSION 

To summarize, our theory describes exactly the thresh- 
old behavior of the CLL correlators. In contrast, the 
'effective exponent theory— strongly distorts the interac- 
tion, making it neither short-ranged nor Coulomb one. 
This is the main reason of the discrepancy between our 
results. 
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APPENDIX A: QUASI-ONE-DIMENSIONAL 
COULOMB INTERACTION POTENTIAL 

Here we derive the interaction potential of electrons, 
confined in a ID quantum wire and interacting via 
Coulomb long-ranged forces. The pair interaction po- 
tential equals 



U(xi,x 2 ) 



' J G(ri - r2)p(r 1 ±)p(r 2 ±)dr 1 ±dr2A 



(Al) 

The a;-axis is directed along the wire, = 
Ti± = (yi,Zi). p{r±) is the transverse density distri- 
bution in the wire, with only first subband populated. 
The transverse distribution is supposed to be localized 
on scale of r± ss d. The Green function is 



G(r) 



1 47T 



(A2) 



Denote x — x\ — xi. The long-distance (x — > oo) be- 
havior of the potential can be well approximated by 
U{x) ps e 2 j\Jx 2 + d 2 . Although this simple approxima- 
tion leads to correct logarithmic behavior of the Fourier- 
transformed potential in the long-wave region of q — > 0, 
it is insufficient to obtain correct expression in the short- 
range region. Below one finds a general derivation. 

Represent the Green function in Eq. (|A1(I via the 
Fourier integral to get 



U(q) 



dxdk 

2^ 2 k 2 ' 



,iqx -ik(ri-r 2 ; 



p(r 1± )p(r2±)dr 1± dr 2 ±. 

(A3) 



Performing integration w.r.t. x gives <5-function, which 
eliminates integration w.r.t. k x . Integration w.r.t. rjj_ 
yields 



U(q) = e 2 



dki 



1 



I" q 2 + k]_ 
For cylindric wire this becomes 

2kdk 



|p(kj 



(A4) 



U(q) = e< 



q 2 + k 2 



\p(k)Y 



(A5) 



Use the fact that p{k) has a scale A ps 1/d, such that 
p(k) is small at k > A. Neglecting for simplicity the 
/c-dependence of p(k) for < k < A, one gets 



U(q) = e 2 



2kdk 

o q 2 + k 2 



e 2 In 



(A6) 



1 

qd 



At qd <ti 1 , this gives the logarithmic divergency 
U(q) fa 2e 2 ln 

When > 1, 

A 2 



(A7) 



(A8) 
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